Giant Wave-Drag Enhancement of Friction in Sliding Carbon Nanotubes 
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Molecular dynamics simulations of coaxial carbon nanotubes in relative sliding motion reveal a 
striking enhancement of friction when phonons whose group velocity is close to the sliding velocity 
of the nanotubes are resonantly excited. The effect is analogous to the dramatic increase in air drag 
experienced by aircraft flying close to the speed of sound, but differs in that it can occur in multiple 
velocity ranges with varying magnitude, depending on the atomic level structures of the nanotubes. 
The phenomenon is a general one that may occur in other nanoscale mechanical systems. 
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The miniaturization of mechanical devices has evolved 
to the point that microelectromechanical systems 
(MEMS) are now in widespread use in a growing num- 
ber of commercial products. Miniaturization has many 
advantages including increased speed and sensitivity, re- 
duced power consumption, and improvements in the po- 
sitional and spatial precision with which complex func- 
tions can be performed. At all length scales the scale- 
relative properties of mechanical systems vary with their 
size due to increasing edge to surface and surface to vol- 
ume ratios. However, as the nanometer length scale is 
approached, an entirely new set of phenomena are be- 
coming relevant that provide great challenges to the fur- 
ther miniaturization of devices. The precise atomic and 
electronic structure of devices, the chemical reactivity of 
device components, the quantum mechanical behaviour 
of atoms, and the presence of large thermal fluctuations 
are all issues of increasing importance at this scale. Engi- 
neering useful nanoscale mechanical devices will require 
an understanding of these effects and their impacts on 
device efficiency and reliability. 

Multiwalled carbon nanotubes (MWCNTs) can be fab- 
ricated with near structural perfection and their resis- 
tance to deformation and the smoothness with which 
coaxial tubes can move relative to one another suggest 
that they could be ideal materials from which to build 
nanoscale machinery. While useful nano-mechanical sys- 
tems have not yet been constructed, some prototype me- 
chanical elements with sub-micron dimensions have been 
i—i 

made using MWCNTs 1 and it seems unlikely that the 
ultimate lower size limit of functional mechanical devices 
has, on the micron scale, already been reached. The- 
ory and simulation can help to investigate the limitations 
of nanoscale machinery and to improve their design by 
providing a more detailed understanding of energy loss 
mechanisms than can be inferred from experiment. This 
Letter is concerned with one of the most important gross 
features of nanoscale phononic friction: Its dependence 
on velocity. 
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FIG. 1: (a) Sliding velocity, v s , as a function of time for 
four different molecular dynamics (MD) simulations of a (5,5) 
CNT sliding within a (10,10) CNT. (b) Average friction force, 
Ff, as a function of the average sliding velocity in near- 
constant-velocity MD simulations. For the (16,2)-(5,5) sys- 
tem, Ff has been scaled by a factor of two for visibility. 



ratios 0, E] mean that components of nanoscale devices 
can move with enormous accelerations and velocities. 
Simulations of an oscillator composed of double-walled 
carbon nanotubes (CNT) |4| showed that even the weak 
van der Waals interactions between layers of a MWCNT 
are sufficient to accelerate nanotubes to velocities of 
thousands of meters per second^- 

In this work, some results of molecular dynamics (MD) 
simulations of CNT "shuttles" , of varying length, slid- 
ing within a larger radius CNT under periodic bound- 
ary conditions Q are reported. The focus is on a striking 
phenomenon, involving an enhancement of friction of sev- 
eral orders of magnitude, that occurs in certain velocity 
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FIG. 2: The difference in (w, fc) space between the power 
spectra of the (10,10) CNT calculated at sliding velocities 
of 1100 ms _1 and 1000 ms _1 . Inset: Frequency-dependent 
power spectrum at different sliding velocities. 
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FIG. 3: The lines are the phonon bands of an isolated (10, 10) 
nanotubes at K. The RB and LA modes are indicated by 
dashed lines. The underlying color plot is the fc-resolved 
phonon power spectrum of the (10, 10) CNT when the (5,5) 
shuttle is moving inside it at 1000 ms" . Circles indicate the 
positions of the peaks seen in Fig|2]while the boxes are points 
that satsify the proposed conditions for enhanced friction. 



ranges and that may have general importance for the de- 
sign of nanomechanical systems 0- Figure ^a) is a plot 
of sliding velocity, v s , as a function of time in MD simula- 
tions of a 8 nm (5, 5) CNT shuttle sliding within a (10, 10) 
CNT in a periodically repeated unit cell of length ~ 22.1 
nm 9]. When v s is in the approximate range 1050 ms" 1 
to 1150ms" 1 , the rate at which the shuttle slows down is 
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FIG. 4: Average over time and azimuthal angle of the local 
radius, r avc of the (10, 10) CNT as a function of displacement 
from one end of the (5, 5) CNT at three different sliding ve- 
locities 



dramatically larger than when v s is outside this range. 

To investigate this effect, MD simulations have been 
performed at close to constant velocity as follows: With 
the CNTs at rest, the system was prepared at a temper- 
ature of 300 K and then equilibrated in a microcanonical 
ensemble for tens of picoseconds. The shuttle was then 
given an initial sliding velocity, v®, along the common 
axis of the nanotubes and MD was performed in a mi- 
crocanonical ensemble, but with v 8 rescaled to v® every 
400 femtoseconds. The average of the net force between 
the two tubes is the average friction force, Ff. When 
the temperature (averaged over 0.5 ps) of either CNT 
rose above the simulation was abandoned and a 

new simulation begun using a different, and independent, 
equilibrated starting configuration. Ff was converged by 
averaging over an aggregate time domain of one or two 
nanoseconds. Extensive tests were performed to verify 
that the results reported here are insensitive to details 
of these simulations such as the frequency of rescaling of 
v s , the tolerance (of 5%) on temperature changes , the 
equilibration time, and the suddenness with which the 
sliding velocity was imposed on the shuttle and subse- 
quently rescaled. 

Figure ^b) is a plot of Ff versus v s . The value of v s 
that is used here is the average velocity in each simu- 
lation, and not v®. As expected from Fig. da), there 
is a massive increase in Ff near v s = 1100 ms" 1 for 
the (10,10)-(5,5) system. Figure QJb) also shows Ff vs 
v s (scaled by a factor of two) for a (5,5) shuttle sliding 
within a (chiral) (16, 2) CNT. In this case the shuttle and 
the outer tube are incommensurate and there are peaks 
in Ff at many different velocities. Although the peaks in 
the Ff are suggestive of resonances, and the possibility 
of resonant friction between CNTs has previously been 
proposed^]] j it wm be shown that this effect can not 
be explained by a resonant mechanism alone. Due to its 
relative simplicity the analysis in the remainder of this 
Letter is focussed on the (10, 10) — (5, 5) system. 

When surfaces slide over one another the relative mo- 
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tion of the corrugated surface potentials directly excites 
particular phonons. A phonon excited in this way can 
be resonantly excited if its frequency coincides with the 
"washboard frequency", w w b = 2nv s /a, where v s is the 
sliding velocity, and a is a length scale that is common 
to both surfaces and that plays a role in the excitation of 
that particular phonon. uj w ^ can take any value depend- 
ing on the value of the v s . 

To identify which phonons are involved in the observed 
enhancement of friction the displacement- dependent 
velocity- velocity correlation function (i.e. A(zi — Zj,t\ — 
£2) = (Vi(ti)vj(t2)), where Zi,Zj are the displacements 
along the axis of the outer nanotube of the primitive unit 
cells to which atoms i and j belong, and Ui(ii), Vjit-z) are 
their velocities at times t\ and t%), of the (10,10) tube 
was Fourier transformed in time and space. This gives 
the frequency (u>) and wavevector (fc = 2-k/X) depen- 
dence of its phonon power spectrum. Figure [5] shows the 
difference in (w, k) space between the power spectra cal- 
culated at 1100 ms _1 and at 1000 ms -1 . One of the two 
peaks in this difference spectrum is centered at (u), A) = 
(146 cm , 00) and the other, less symmetric peak, has 
maximum intensity at (u>, A) = (140 cm -1 , 3.68 nm) and 
a substantial shoulder at (137 cm -1 , 3.16nm). The in- 
set to Fig|21 demonstrates that, apart from these modes 
at u> ~ 140 cm" 1 that are strongly excited at v s — 
1100 ms _1 , the phonon power spectrum is rather insen- 
sitive to the sliding velocity. 

The phonon band structure of the (10, 10) CNT at 
K was calculated by diagonalizing its dynamical matrix. 
These bands are superimposed in Fig. [3] on an image 
depicting the fc-resolved power spectrum calculated from 
the MD simulations at 1000 ms" 1 . The phonon energet- 
ics are obviously different in the 300 K system with the 
shuttle present, however, some important conclusions can 
be drawn from this data. First of all, there is a wide range 
of velocities at which the washboard frequency coincides 
with phonon frequencies. Even if attention is restricted 
to the radial breathing (RB) mode, "washboard reso- 
nances" can occur at many velocities between 1000 ms _1 
and 2500 ins -1 . This means that resonance alone cannot 
account for the observed dependence of Ff on v s . 

At this point we postulate that increased friction oc- 
curs when there is resonant excitation of a phonon in the 
outer nanotube that moves at close to the velocity of the 
shuttle. In other words, the group velocity, v g , of the 
phonon that is excited by the sliding motion is close to 
the value of v s required for a match between the phonon 
frequency, u>(k), and w w b, he. the following two crite- 
ria are simultaneously met: (i) v g (k) = duj(k)/dk w v s 
and, (ii)w w b(v) ~ ui(k) =>■ w w b = w w ba/27r rj v s . Evi- 
dence to support this hypothesis will now be provided. 
This evidence consists of a demonstration that these cri- 
teria can only be met at velocities close to 1100 ms" 1 by 
phonons that have frequencies and wavevectors consis- 
tent with the anomalous phonon excitation demonstrated 



in FigEl an d furthermore, that these phonon modes are 
those that would, most obviously, be excited by the wash- 
board mechanism - the axially symmetric radial breath- 
ing (RB) and longitudinal acoustic (LA) modes. The LA 
mode and the RB mode both have perfect axial symme- 
try and they mix strongly with each other and with other 
phonon modes near (oj, k) = (140, 0.1) where a band anti- 
crossing occurs. 

In FigI3 circles mark the positions in (lj, k) space of 
the peaks seen in Fig|3 Four other points are marked 
with boxes. At these points criteria (i) and (ii) are both 
met. Criteria (i) and (ii) are also met at a number of 
other points |Tll|. however, none of these other phonons 
either have axial symmetry or group velocities in the 
range 950ms _1 < v g < 1250ms -1 . Furthermore, of these 
phonons, all but one (for which u> rs 128cm _1 ) have fre- 
quencies outside the range 120cm _1 < uj < 160cm -1 . 
The four marked points in Fig. [3] mostly have the char- 
acter of RB modes with some LA character in the points 
closest to the band anti-crossing. 

The agreement between two of these four phonons and 
the simulation results presented in figures ^ an d [21 is re ~ 
markably good given that they have been calculated from 
the K phonons of an isolated (10, 10) tube. The val- 
ues of v s and i; w b are close to 1100 ms _1 and they are 
close in (u>, k) space to the positions at which anomalous 
phonon excitation occurs in the MD simulations. Two 
other phonons do not, at first, appear consistent with 
figures H and The first of these has a group velocity 
of ~ 500 ms _1 and no peak is observed at this veloc- 
ity in Fig. ^ However, the wavelength of this phonon 
is A w 24.5 nm which is larger than the length of the 
simulation cell, and therefore it cannot be excited in the 
MD simulations. The second point, at A « 1.5 nm, has 
a group velocity close to 1100 ms -1 but no anomalous 
phonon excitation at A « 1.5 nm is observed in the MD 
simulations. A possible reason for this is provided below. 

The MD simulations have been repeated with shut- 
tles of lengths (in nm, and excluding the lengths of the 
caps) 10, 8, 5.3, 2.6 and (i.e. a Ceo molecule). No en- 
hancement of Ff was observed near 1100 ms -1 or at any 
other velocity for the Ceo molecule. There were peaks 
at u s = 1100 ms -1 for all other shuttle lengths and 
the heights of the peaks increased almost linearly with 
the shuttle length. For the 2.6 nm shuttle, there was 
anomalous phonon excitation at the smaller wavelength 
A Ri 3.2 — 3.7 nm but no anomalous phonon excitation 
at A — > 00. As the length of the shuttle was increased, 
the A — > 00 peak began to grow at the expense of the 
smaller wavelength peak until, for the longest shuttle, 
this smaller wavelength peak had almost disappeared. It 
is reasonable that short shuttles cannot excite very long 
wavelength phonons and that, as the shuttle length in- 
creases, phonons with a larger wavelength are excited. 
Furthermore, when the washboard excitation is in-phase 
with a phonon at a given point along the outer tube 
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it must be in anti-phase with a point a distance A/2 
away. Therefore, a shuttle cannot resonantly excite a 
phonon whose wavelength is substantially smaller than 
twice its length. This explains the shift in intensity of 
the anomalous phonon excitation from shorter to longer 
wavelengths as the shuttle length increases, and the ab- 
sence of anomalous excitation of phonons at a wavelength 
of ~ 1.5 nm. 

A physical mechanism is now proposed for the results 
of the MD simulations: When v s is sufficiently close to 
the group velocity of a RB mode of the outer nanotube 
that is being resonantly excited by the relative motion 
of the corrugated surface potentials of the nanotubes, a 
wave packet remains in the vicinity of the shuttle for long 
enough that it is reinforced and amplified. The build up 
of a radial distortion in the outer nanotube that moves 
along with the shuttle results. This distortion, in turn, 
changes the interaction between the two nanotubes - a 
feedback that makes the response of the outer nanotube 
to the presence of the shuttle strongly non-linear. The 
details of this non-linearity are likely to be complicated 
and sensitively velocity-dependent, but the general effect 
is analogous to the sharp increase in drag experienced by 
airplanes travelling at close to the speed of sound. At low 
speeds, an airplane creates a disturbance in the air and 
sends out pressure pulses that travel ahead of it to sepa- 
rate the air flow, thereby allowing air to move smoothly 
around it. As Mach 1 is approached, however, pressure 
waves build up to form a shock front just ahead of the 
airplane. Air in the airplane's path has little warning of 
its arrival and so the air's compressibility plays a cru- 
cial role in the flow dynamics. The dynamics become 
strongly nonlinear, instabilities appear in the flow, and 
lots of energy is dissipated by the shock front. The drag 
that results is commonly known as "wave drag" ^ an( i 
its maximum as a function of velocity is the "sound bar- 
rier" . At supersonic speeds the flow stabilizes and the 
drag is reduced. 

Further evidence that a similar mechanism is respon- 
sible for the enhancement of friction between CNTs is 
provided in Fig0] which shows the average radius, r avo , 
of the outer CNT as a function of displacement, Z, from 
one end of the shuttle. The average is performed over 
the azimuthal angle and over time. At velocities of 
500 ms -1 and 1000 ms -1 , r avc is very similar, however, 
at 1100 ms" 1 there is a large change, suggesting that a 
wave builds up around the shuttle at this velocity. The 
fact that the difference between r ave at 1100 ms -1 and at 
1000 ms -1 is of the same order of magnitude as the vari- 
ations in r avc at either 500 ms" 1 or 1000 ms -1 implies 
that the inter-tube interaction is strongly affected by the 
distortion. This, in turn, implies a non-linear response 
of the outer CNT to the shuttle that could explain the 
sharp increase in the rate of energy dissipation. 

A vast literature exists on wave-drag and related phe- 
nomena in the context of aeronautics 121 and much of the 



theory may be adaptable to the nanoscale. One impor- 
tant difference at the nanoscale, however, is that, depend- 
ing on the atomic level structure, this effect can occur in 
the same device at multiple velocities. Particularly for 
incommensurate systems, there may be more than one 
length scale that is common to the contacting surfaces 
and that can cause washboard resonances. Furthermore, 
multiple phonons with different group velocities may sat- 
isfy the criteria necessary for this effect to occur. The 
(16, 2) — (5, 5) system shown in figure^b) is an example 
in which there are many peaks in Ff vs v s . 

In conclusion, a mechanism by which huge increases in 
mechanical energy dissipation at well defined velocities 
can occur at the nanoscale has been demonstrated using 
atomistic simulations of coaxial carbon nanotubes sliding 
relative to one another. As demonstrated in Fig. Qfa), 
this effect can have a dramatic impact on the dynamics 
of a nano-mechanical system. 
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